Aldehyde dehydrogenase diversity in bacteria of the Pseudomonas genus.

Aldehyde dehydrogenases (ALDHs) comprise one of the most ancient protein superfamilies widely distributed in the three domains of life. Their members have been extensively studied in animals and plants, sorted out in different ALDH protein families and their participation in a broad variety of metabolic pathways has been documented. Paradoxically, no systematic studies comprising ALDHs from bacteria have been performed so far. Among bacteria, the genus Pseudomonas occupies numerous ecological niches, and is one of the most complex bacterial genera with the largest number of known species. For these reasons, we selected Pseudomonas as a paradigm to analyze the diversity of ALDHs in bacteria. With this aim, complete Pseudomonas genome sequences and annotations were retrieved from NCBI's RefSeq genome database. The 258 Pseudomonas strains belong to 46 different species, along with 23 with no species designation. The genomes of these Pseudomonas contain from 3,315 to 6,825 annotated protein coding genes. A total of 6,510 ALDH sequences were found in the selected Pseudomonas, with a median of 24 ALDH-coding genes per strain (by comparison humans possess only 19 different ALDH loci). Pseudomonas saudiphocaensis possesses the lowest number of aldh genes (9), while Pseudomonas pseudoalcaligenes KF707 NBRC110670 possesses the maximum number of aldh genes (49). The ALDHs found in Pseudomonas can be sorted out into 42 protein families, with a predominance of 14 families, which contained 76% of all ALDHs found. In this regard, it is important to note that many Pseudomonas genomes have multiple aldh genes coding for proteins belonging to the same family. Given that all strains contained members of families ALDH4, ALDH5, ALDH6, ALDH14, ALDH18 and ALDH27, we consider these families to be part of the core Pseudomonas genome.

Correction: Diversity and Evolutionary Analysis of Iron-Containing (Type-III) Alcohol Dehydrogenases in Eukaryotes.

[This corrects the article DOI: 10.1371/journal.pone.0166851.].

Diversity and Evolutionary Analysis of Iron-Containing (Type-III) Alcohol Dehydrogenases in Eukaryotes.

BACKGROUND: Alcohol dehydrogenase (ADH) activity is widely distributed in the three domains of life. Currently, there are three non-homologous NAD(P)+-dependent ADH families reported: Type I ADH comprises Zn-dependent ADHs; type II ADH comprises short-chain ADHs described first in Drosophila; and, type III ADH comprises iron-containing ADHs (FeADHs). These three families arose independently throughout evolution and possess different structures and mechanisms of reaction. While types I and II ADHs have been extensively studied, analyses about the evolution and diversity of (type III) FeADHs have not been published yet. Therefore in this work, a phylogenetic analysis of FeADHs was performed to get insights into the evolution of this protein family, as well as explore the diversity of FeADHs in eukaryotes. PRINCIPAL FINDINGS: Results showed that FeADHs from eukaryotes are distributed in thirteen protein subfamilies, eight of them possessing protein sequences distributed in the three domains of life. Interestingly, none of these protein subfamilies possess protein sequences found simultaneously in animals, plants and fungi. Many FeADHs are activated by or contain Fe2+, but many others bind to a variety of metals, or even lack of metal cofactor. Animal FeADHs are found in just one protein subfamily, the hydroxyacid-oxoacid transhydrogenase (HOT) subfamily, which includes protein sequences widely distributed in fungi, but not in plants), and in several taxa from lower eukaryotes, bacteria and archaea. Fungi FeADHs are found mainly in two subfamilies: HOT and maleylacetate reductase (MAR), but some can be found also in other three different protein subfamilies. Plant FeADHs are found only in chlorophyta but not in higher plants, and are distributed in three different protein subfamilies. CONCLUSIONS/SIGNIFICANCE: FeADHs are a diverse and ancient protein family that shares a common 3D scaffold with a patchy distribution in eukaryotes. The majority of sequenced FeADHs from eukaryotes are distributed in just two subfamilies, HOT and MAR (found mainly in animals and fungi). These two subfamilies comprise almost 85% of all sequenced FeADHs in eukaryotes.

Residues that influence coenzyme preference in the aldehyde dehydrogenases.

To find out the residues that influence the coenzyme preference of aldehyde dehydrogenases (ALDHs), we reviewed, analyzed and correlated data from their known crystal structures and amino-acid sequences with their published kinetic parameters for NAD(P)(+). We found that the conformation of the Rossmann-fold loops participating in binding the adenosine ribose is very conserved among ALDHs, so that coenzyme specificity is mainly determined by the nature of the residue at position 195 (human ALDH2 numbering). Enzymes with glutamate or proline at 195 prefer NAD(+) because the side-chains of these residues electrostatically and/or sterically repel the 2'-phosphate group of NADP(+). But contrary to the conformational rigidity of proline, the conformational flexibility of glutamate may allow NADP(+)-binding in some enzymes by moving the carboxyl group away from the 2'-phosphate group, which is possible if a small neutral residue is located at position 224, and favored if the residue at position 53 interacts with Glu195 in a NADP(+)-compatible conformation. Of the residues found at position 195, only glutamate interacts with the NAD(+)-adenosine ribose; glutamine and histidine cannot since their side-chain points are opposite to the ribose, probably because the absence of the electrostatic attraction by the conserved nearby Lys192, or its electrostatic repulsion, respectively. The shorter side-chains of other residues-aspartate, serine, threonine, alanine, valine, leucine, or isoleucine-are distant from the ribose but leave room for binding the 2'-phosphate group. Generally, enzymes having a residue different from Glu bind NAD(+) with less affinity, but they can also bind NADP(+) even sometimes with higher affinity than NAD(+), as do enzymes containing Thr/Ser/Gln195. Coenzyme preference is a variable feature within many ALDH families, consistent with being mainly dependent on a single residue that apparently has no other structural or functional roles, and therefore can easily be changed through evolution and selected in response to physiological needs.

Amino acid residues that affect the basicity of the catalytic glutamate of the hydrolytic aldehyde dehydrogenases.

In the catalytic mechanism of hydrolytic aldehyde dehydrogenases (ALDHs) the role of Glu268 (mature human ALDH2 numbering) as a general base is of major relevance. Since Glu268 basicity depends on its protein environment, here we explore its interactions with other amino acid residues in the three different conformations observed in ALDH crystal-structures: "inside", "intermediate" and "outside". In all of them Glu268 is in a hydrophobic environment. In the "inside" conformation, the theoretical pKa estimated by PROPKA3 is the result of the effects of hydrogen bonds with the protonated thiol of the catalytic Cys302 and/or the main-chain amide nitrogen of the highly conserved Gly270, and of charge-charge interactions with neighboring side-chains-Lys178, Glu/Asp476, His465 or Glu399 depending on the enzyme. In the "intermediate" conformation Glu268-carboxyl pKa is influenced by interactions with Glu/Asp476, Arg/Lys475, Lys/Arg178, His465 or Arg459, also depending on the enzyme. In the "outside" conformation, the effects on Glu268-carboxyl pKa arise from hydrogen bonds with the side chains of the strictly conserved Thr224 and/or of Lys/Arg178, and from charge-charge interactions with Lys/Arg/Asp178, Glu476, or Arg459. The estimated pKas and interactions of Glu268-carboxyl in the "intermediate" and "outside" conformations are consistent with their previously proposed roles in activating the hydrolytic water and in a proton relay mechanism, respectively. Water channels connecting Glu268 with the bulk water were found in all hydrolytic ALDHs. In the "inside" conformation the theoretical pKas of the Glu268-carboxyl and Cys302-thiol groups suggest that the carboxyl cannot receive the proton from the thiol. We propose that a protonated Cys302 might perform the nucleophilic attack on the aldehyde, which can be facilitated by Glu268 in the "intermediate" conformation. Finally, the conservation of the residues influencing Glu268 basicity between and within ALDH families suggests that these residues, not previously studied, are important for the catalytic mechanism of many ALDH enzymes.

Exploring the evolutionary route of the acquisition of betaine aldehyde dehydrogenase activity by plant ALDH10 enzymes: implications for the synthesis of the osmoprotectant glycine betaine.

BACKGROUND: Plant ALDH10 enzymes are aminoaldehyde dehydrogenases (AMADHs) that oxidize different ω-amino or trimethylammonium aldehydes, but only some of them have betaine aldehyde dehydrogenase (BADH) activity and produce the osmoprotectant glycine betaine (GB). The latter enzymes possess alanine or cysteine at position 441 (numbering of the spinach enzyme, SoBADH), while those ALDH10s that cannot oxidize betaine aldehyde (BAL) have isoleucine at this position. Only the plants that contain A441- or C441-type ALDH10 isoenzymes accumulate GB in response to osmotic stress. In this work we explored the evolutionary history of the acquisition of BAL specificity by plant ALDH10s. RESULTS: We performed extensive phylogenetic analyses and constructed and characterized, kinetically and structurally, four SoBADH variants that simulate the parsimonious intermediates in the evolutionary pathway from I441-type to A441- or C441-type enzymes. All mutants had a correct folding, average thermal stabilities and similar activity with aminopropionaldehyde, but whereas A441S and A441T exhibited significant activity with BAL, A441V and A441F did not. The kinetics of the mutants were consistent with their predicted structural features obtained by modeling, and confirmed the importance of position 441 for BAL specificity. The acquisition of BADH activity could have happened through any of these intermediates without detriment of the original function or protein stability. Phylogenetic studies showed that this event occurred independently several times during angiosperms evolution when an ALDH10 gene duplicate changed the critical Ile residue for Ala or Cys in two consecutive single mutations. ALDH10 isoenzymes frequently group in two clades within a plant family: one includes peroxisomal I441-type, the other peroxisomal and non-peroxisomal I441-, A441- or C441-type. Interestingly, high GB-accumulators plants have non-peroxisomal A441- or C441-type isoenzymes, while low-GB accumulators have the peroxisomal C441-type, suggesting some limitations in the peroxisomal GB synthesis. CONCLUSION: Our findings shed light on the evolution of the synthesis of GB in plants, a metabolic trait of most ecological and physiological relevance for their tolerance to drought, hypersaline soils and cold. Together, our results are consistent with smooth evolutionary pathways for the acquisition of the BADH function from ancestral I441-type AMADHs, thus explaining the relatively high occurrence of this event.

An Lrp-type transcriptional regulator controls expression of the Bacillus subtilis chromate transporter.

The Bacillus subtilis strain 168 genome contains the chr3N-chr3C genes encoding the Chr3N/Chr3C protein pair of the chromate ion transporter (CHR) superfamily. Chr3N/Chr3C confers chromate resistance in Escherichia coli only when both proteins are expressed. Upstream of chr3N is the chrS gene encoding ChrS, a protein with homology to the Lrp/AsnC family of transcriptional regulators. When the chrS-chr3N-chr3C gene cluster was transferred to E. coli, a diminished level of chromate resistance was observed, as compared with E. coli transformants bearing only the chromate resistance genes, which displayed full resistance. These data suggested that the chrS gene product acts as negative regulator. RT-PCR assays demonstrated that expression of chrS diminishes transcription of the chromate resistance genes in E. coli, and that this repression was overcome by chromate. Electrophoretic mobility shift assays showed that purified ChrS protein specifically binds to the 5' region of chrS. These results indicate that the chr gene cluster forms an operon regulated negatively by ChrS binding to its own gene's regulatory region, and positively by chromate ions. Sequence analysis revealed similar operons in many Bacillales strains, suggesting some adaptive advantage. This is the first example of a bacterial heavy-metal resistance system controlled by an Lrp-type transcriptional regulator.

Potential monovalent cation-binding sites in aldehyde dehydrogenases.

Potassium ions are non-essential activators of several aldehyde dehydrogenases (ALDHs), whereas a few others require the cation for activity. Two kinds of cation-binding sites, which we named intra-subunit and inter-subunit, have been observed in crystal structures of ALDHs, and based on reported crystallographic data, we here propose the existence of a third kind located in the central cavity of some tetrameric ALDHs. Given the high structural similarity between these enzymes, cation-binding sites may be present in many other members of this superfamily. To explore the prevalence of these sites, we compared 37 known crystal structures from 13 different ALDH families and evaluated the possible existence of a cation on the basis of the number, distance and geometry of its potential interactions, as well as of B-factor values of modeled cations obtained in new refinements of some reported crystal structures. Also, by performing multiple alignments of 855 non-redundant amino acid sequences, we assessed the degree of conservation in their respective families of the amino acid residues putatively relevant for cation binding. Among the ALDH enzymes studied, and according to our analyses, potential intra-subunit cation-binding sites seem to be present in most members of ALDH2, ALDH1L, ALDH4, ALDH5, ALDH7, ALDH10, and ALDH25 families, as well as in the bacterial and fungal members of the ALDH9 family and in a few ALDH1, ALDH6, ALDH11 and ALDH26 enzymes; potential inter-subunit sites in members of ALDH1L, ALDH3, ALDH4 from bacillales, ALDH5, ALDH7, ALDH9, ALDH10, ALDH11 and ALDH25 families; and potential central-cavity sites only in some bacterial and animal ALDH9s and in most members of the ALDH1L family. Because potassium is the most abundant intracellular cation, we propose that these are potassium-binding sites, but the specific structural and/or functional roles of the cation bound to these different sites remain to be investigated.

Structural determinants of substrate specificity in aldehyde dehydrogenases.

Within the aldehyde dehydrogenase (ALDH) superfamily, proteins belonging to the ALDH9, ALDH10, ALDH25, ALDH26 and ALDH27 families display activity as ω-aminoaldehyde dehydrogenases (AMADHs). These enzymes participate in polyamine, choline and arginine catabolism, as well as in synthesis of several osmoprotectants and carnitine. Active site aromatic and acidic residues are involved in binding the ω-aminoaldehydes in plant ALDH10 enzymes. In order to ascertain the degree of conservation of these residues among AMADHs and to evaluate their possible relevance in determining the aminoaldehyde specificity, we compared the known amino acid sequences of every ALDH family that have at least one member with known crystal structure, as well as the electrostatic potential surface of the aldehyde binding sites of these structures. Our analyses showed that four or three aromatic residues form a similar "aromatic box" in the active site of the AMADH enzymes, being the equivalents to Phe170 and Trp177 (human ALDH2 numbering) strictly conserved in all of them, which supports their relevance in binding the aminoaldehyde by cation-π interactions. In addition, all AMADHs exhibit a negative electrostatic potential surface in the aldehyde-entrance tunnel, due to side-chain carboxyl and hydroxyl groups or main-chain carbonyl groups. In contrast, ALDHs that have non-polar or negatively charged substrates exhibit neutral or positive electrostatic potential surfaces, respectively. Finally, our comparative sequence analyses revealed that the residues equivalent to Asp121 and Phe170 are highly conserved in many ALDH families irrespective of their substrate specificity-suggesting that they perform a role in catalysis additional or different to binding of the substrate-and that the positions Met124, Cys301, and Cys303 are hot spots changed during evolution to confer aldehyde specificity to several ALDH families.

Natural alcohol exposure: is ethanol the main substrate for alcohol dehydrogenases in animals?

Alcohol dehydrogenase (ADH) activity is widely distributed in all phyla. In animals, three non-homologous NAD(P)(+)-dependent ADH protein families are reported. These arose independently throughout evolution and possess different structures and mechanisms of reaction: type I (medium-chain) ADHs are zinc-containing enzymes and comprise the most studied group in vertebrates; type II (short-chain) ADHs lack metal cofactor and have been extensively studied in Drosophila; and type III ADHs are iron-dependent/-activated enzymes that were initially identified only in microorganisms. The presence of these different ADHs in animals has been assumed to be a consequence of chronic exposure to ethanol. By far the most common natural source of ethanol is fermentation of fruit sugars by yeast, and available data support that this fruit trait evolved in concert with the characteristics of their frugivorous seed dispersers. Therefore, if the presence of ADHs in animals evolved as an adaptive response to dietary ethanol exposure, then it can be expected that the enzymogenesis of these enzymes began after the appearance of angiosperms with fleshy fruits, because substrate availability must precede enzyme selection. In this work, available evidence supporting this possibility is discussed. Phylogenetic analyses reveal that type II ADHs suffered several duplications, all of these restricted to flies (order Diptera). Induction of type II Adh by ethanol exposure, a positive correlation between ADH activity and ethanol resistance, and the fact that flies and type II Adh diversification occurred in concert with angiosperm diversification, strongly suggest that type II ADHs were recruited to allow larval flies to exploit new restricted niches with high ethanol content. In contrast, phyletic distribution of types I and III ADHs in animals showed that these appeared before angiosperms and land plants, independently of ethanol availability. Because these enzymes are not induced by ethanol exposure and possess a high affinity and/or catalytic efficiency for non-ethanol endogenous substrates, it can be concluded that the participation of types I and III ADHs in ethanol metabolism can be considered as incidental, and not adaptive.

Phylogenetic analysis of the chromate ion transporter (CHR) superfamily.

ChrA is a membrane protein that confers resistance to the toxic ion chromate through the energy-dependent chromate efflux from the cytoplasm. In the protein databases, ChrA is a member of the chromate ion transporter (CHR) superfamily, composed of at least several dozens of members, distributed in the three domains of life. The aim of this work was to perform a phylogenetic analysis of the CHR superfamily. An exhaustive search for ChrA homologous proteins was carried out at the National Center for Biotechnology Information database. One hundred and thirty-five sequences were identified as members of the CHR superfamily [77 long-chain sequences, or bidomains (LCHR), and 58 short-chain sequences, or monodomains (SCHR)], organized mainly as tandem pairs of genes whose resultant proteins probably possess oppositely oriented membrane topology. LCHR sequences were split into amino and carboxyl domains, and the resultant domains were aligned with the SCHR proteins. A phylogenetic tree was reconstructed using four different methods, obtaining similar results. The domains were grouped into three clusters: the SCHR proteins cluster, the amino domain cluster of LCHR proteins and the carboxyl domain cluster of LCHR proteins. These results, as well as differences in the genomic context of CHR proteins, enabled the proteins to be sorted into two families (SCHR and LCHR), and 10 subfamilies. Evidence was found suggesting an ancient origin of LCHR proteins from the fusion of two SCHR protein-encoding genes; however, some secondary events of fusion and fission may have occurred later. The separate distribution of the LCHR and SCHR proteins, differences in the genomic context in both groups and the fact that chromate transport has been demonstrated only in LCHR proteins suggest that the CHR proteins comprise two or more paralogous groups in the CHR superfamily.

Diversity, taxonomy and evolution of medium-chain dehydrogenase/reductase superfamily.

A comprehensive, structural and functional, in silico analysis of the medium-chain dehydrogenase/reductase (MDR) superfamily, including 583 proteins, was carried out by use of extensive database mining and the blastp program in an iterative manner to identify all known members of the superfamily. Based on phylogenetic, sequence, and functional similarities, the protein members of the MDR superfamily were classified into three different taxonomic categories: (a) subfamilies, consisting of a closed group containing a set of ideally orthologous proteins that perform the same function; (b) families, each comprising a cluster of monophyletic subfamilies that possess significant sequence identity among them and might share or not common substrates or mechanisms of reaction; and (c) macrofamilies, each comprising a cluster of monophyletic protein families with protein members from the three domains of life, which includes at least one subfamily member that displays activity related to a very ancient metabolic pathway. In this context, a superfamily is a group of homologous protein families (and/or macrofamilies) with monophyletic origin that shares at least a barely detectable sequence similarity, but showing the same 3D fold. The MDR superfamily encloses three macrofamilies, with eight families and 49 subfamilies. These subfamilies exhibit great functional diversity including noncatalytic members with different subcellular, phylogenetic, and species distributions. This results from constant enzymogenesis and proteinogenesis within each kingdom, and highlights the huge plasticity that MDR superfamily members possess. Thus, through evolution a great number of taxa-specific new functions were acquired by MDRs. The generation of new functions fulfilled by proteins, can be considered as the essence of protein evolution. The mechanisms of protein evolution inside MDR are not constrained to conserve substrate specificity and/or chemistry of catalysis. In consequence, MDR functional diversity is more complex than sequence diversity. MDR is a very ancient protein superfamily that existed in the last universal common ancestor. It had at least two (and probably three) different ancestral activities related to formaldehyde metabolism and alcoholic fermentation. Eukaryotic members of this superfamily are more related to bacterial than to archaeal members; horizontal gene transfer among the domains of life appears to be a rare event in modern organisms.

Ozone and Its Precursors in the Atmosphere of Mexico City.

The relative apportionment of hydrocarbons (HCs) coming from mobile, fixed, and other sources (not correlated either to carbon monoxide [CO] or sulfur dioxide [SO2] emissions as solvent evaporation and biogenic sources) is calculated as previously proposed by Riveros et al.1 through the linear approximation [HC]tol = [HC]0 + m1 [CO] + m2 [SO2], where m1 and m2 are fitted constants. The obtained apportionment with air samples taken in 1993 is consistent with the earlier published apportionment with air samples taken in 1992, validating the previous procedure. This analysis suggests that 75% of HC originate from mobile sources, 5-18% from fixed sources, and 7-20% from other sources (mainly solvents and bio-genic sources). A similar analysis was employed to estimate the relative contribution of HCs and nitric oxides (NO2) to ozone (O3) formation, the most important air pollutant in Mexico City. In this way, through a local lineation of O3 isopleths, simultaneous measurements of HC and NO2 in the atmosphere were fitted to the equation-[O3]peak = [O3]0 + ma [HC] + mb [NO2]-to predict O3 peak. With these data, the adjusted parameters show that NO2, not HC as was proposed previously, is the most important contributor to O3 formation.

Enzymology of ethanol and acetaldehyde metabolism in mammals

Ehtanol or wthyl alcohol is a molecule that, in mammals, is naturally present at low concentrations due to its production by gastrointestinal flora fermentation activity. However, it is remarkable that this metabolite, with a clearly minor role in regular vertebrate metabolism, can be oxidized into acetaldehyde through several ensymatic and non-enzymatic mechanisms, which comprise the activity of more than ten ensymes and isozymes, many of them broadly distributed in different specie and tissues. In correspondence, acetaldehyde can also be oxidized into acetate through several enzymatic pathways that involve about ten enzymes and isozymes which also have a broad distribution In this article, a complete review of the aforementioned metabolic pathways is elaborated. From this group, the participation and wide distribution of alcohol dehydrogenase and aldehyde dehydrogenase systems are emphasized. The mechanism of reaction, kinetic characteristics and physiological relevance are described, and finally, the possible physiological role of these enzymatic systems as responsible to synthesize or catabolize several endogenous metabolites that regulate growth, metabolism, differentiation and neuroendocrine function in mammals are discussed